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\f} • Abstract 
t— i ; 

We consider stability under evolution as a criterion to select a physical boundary 
q ' state for the spinfoam formalism. As an example, we apply it to the simplest spin- 

foam defined by a single quantum tetrahedron and solve the associated eigenvalue 
5— i [ problem at leading order in the large spin limit. We show that this fixes uniquely 

the free parameters entering the boundary state. Remarkably, the state obtained 
this way gives a correlation between edges which runs at leading order with the 
inverse distance between the edges, in agreement with the linearized continuum 
' theory. Finally, we give an argument why this correlator represents the propagation 

of a pure gauge, consistently with the absence of physical degrees of freedom in 3d 
, general relativity. 

i> : 1 Introduction 

o ■ 

LQG is a background independent approach to the quantization of General Relativity 
^ (GR) [1]. In its covariant spinfoam formulation, it provides us with a regularized ex- 
pression for the path integral of the full theory. The state of the art in spinfoams is 
the proposal for computing the graviton propagator in non-perturbative quantum gravity 
[2, 3, 4, 5, 6, 7, 8]. One works with a bounded region of space-time and introduces a semi- 
classical state peaked on a given boundary geometry. This boundary state gauge-fixes 
the spinfoam amplitude [9] and allows to compute the (two-point) correlations for the 
gravitational field [2], thus effectively inducing a non-trivial bulk geometry. This can be 
seen as the proper way to introduce a background metric in this context of background 
independent quantum gravity. It further allows a perturbative expansion (in the inverse of 
the observation length scale) with a leading order corresponding to a flat spacetime metric 
at large scale. The hope is that this will allow in the long run to perform calculations in 
quantum gravity and extract useful phenomenological predictions. 
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This program has produced a number of interesting results. First, the graviton prop- 
agator point of view allowed a criticism and a detailed description of the shortcomings 
of the Barrett-Crane model [10], which has been the most studied spinfoam model up to 
now [4, 7, 11], and spurred the developments of promising alternatives [12]. Beside this, 
a detailed analysis of the perturbative expansion, with both analytical and numerical 
computations, was performed in the context of 3d quantum gravity [3, 5]. Also many 
mathematical and technical tools have been developed to compute and study the 4d spin- 
foam correlations [4, 6] and there has been a preliminary computation of the three-point 
correlation functions [13]. Finally, there has been a thorough numerical study of the 4d 
correlations for the Barrett-Crane model, which confirmed the semiclassical behaviour 
at large scale but also allowed to study the short scale behaviour and show that the 
correlations get dynamically regularized at the Planck scale [8]. 

The boundary state is a key ingredient of this approach. The original ansatz [2] is a 
Gaussian state with a phase factor in the Hilbert space of boundary spin networks. This 
ansatz has been improved in [5, 6] to allow an easier study of the perturbative expansion. 
These choices do lead to the right behavior for the correlations in the semiclassical limit, 
but leave two important questions open: 

(i) There are free parameters in the Gaussian ansatz. Can we fix these parameters 
dynamically? 

(ii) The Gaussian state can only be a first order approximation of some physical state 
in the full theory. How can we characterize the full boundary state? 

We believe that both questions can be addressed taking seriously the fact that we 
need to work with a physical state, that is gauge invariant spin network states that solve 
the quantum gravity constraints (scalar and momentum constraints). This issue has not 
been addressed yet in the literature. We expect that such a dynamical selection of the 
boundary state will help characterizing it in the full theory, thus addressing point (ii), but 
more importantly we expect it to uniquely fix the parameters of the Gaussian ansatz at 
leading order, thus addressing point (i) and increasing the predicting power of the theory. 

In this paper, we address these issues in 3d Riemannian quantum gravity. The theory 
is indeed much simpler in three spacetime dimensions. Gravity is topological and we know 
how to spinfoam quantize it exactly as the Ponzano-Regge model [14, 15]. We can solve 
this toy model explicitly and we know the physical states (spin network functionals on the 
space of flat connections up to diff eomorphisms) . It thus seems to be the perfect arena 
to address the problem before tackling it in the framework of 4d spinfoam models. 

To keep things simple and explicit, we consider the smallest 3d triangulation, that is 
a single tetrahedron. The corresponding Ponzano-Regge spinfoam amplitude is simply 
Wigner's {6j} symbol for the unitary group SU(2). Following the setting introduced 
in [3, 5], we consider the "time-fixed" tetrahedron: out of its six edges, we study the 
correlations between the fluctuations of two opposite edges while freezing the lengths of 
the remaining four edges. The fixed-length edges define the time interval associated to 
that piece of 3d spacetime, while the two fluctuating edges are distinguished as the initial 
and final edges. Carefully defining a physical state as a wave function unaffected by 
time evolution, we show that the phased Gaussian ansatz actually turns out to be a true 
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physical state at first order in the asymptotical regime. 

More precisely, the physical state criterion uniquely fixes the width of the Gaussian. 
First it fixes the (leading order) scaling of the Gaussian width with the observation length 
scale, which turns out to be the right scaling in order to recover a flat metric in the semi- 
classical regime [3, 5, 16]. Second it fixes the exact factor in front of the scaling, which 
turns out to be such that the leading order of the graviton propagator goes precisely as 
the inverse distance between the edges. 

Finally, we look at the tetrahedron from a covariant point of view. There is a unique 
physical state, namely the flat connection boundary state. We discuss the gauge-fixing 
procedure to go down to the "time-fixed" tetrahedron. We also address the issue that 
there is actually no graviton in 3d gravity (since it is a topological theory without local 
degree of freedom) and argue that the 3d spinfoam graviton correlations entirely come 
from the choice of gauge fixing. 



2 From the propagation kernel to the boundary state 

To better understand the physics behind this paper, let us briefly recall some basic facts 
of quantum mechanics and quantum field theory, concerning the propagation kernel and 
the vacuum state, in addition to their physical interpretation. These objects incorporate 
key aspects of the interpretative problems of background independent quantum gravity. 

Consider a scalar in the Schrodinger picture [17, 18, 19], and introduce two spacelike 
hyperplanes in Minkowski spacetime, separed by a time T. Denote 0i and 02 two classical 
field configurations associated with these planes, and ip n [<p} a complete basis of energy 
eigenstates. Given a state i/)[<f>i] = ^n^^nf^i] on ^ ne initial plane, the evolution to the 
final state can be written in terms of the propagation kernel K[(f>i,(f>2,T\ through 

V# 2 ,T)= f P0 1 ir[0 1 ,0 2 ,T]^(0i,O) = ^c n e-^ T ^ n [0 2 ,O]. (1) 

J n 

The second equality shows that the kernel (that can be evaluated as a path integral with 
0i and 02 as boundary data) can be decomposed as 

x[0i,0 2 ,T] = Y, e ~ iEnT M<i>iWM- (2) 

n 

Among the n-particle states in (2), the vacuum ip [0] is characterized by having the 
minimal energy. It can be extracted directly from the kernel using a Wick rotation to 
Euclidean time T E = iT, 

lim K[<f>,0,-iT E ] =Vo[0]- (3) 

Tei— >oo 

Finally, using the kernel (2) and vacuum (3) one can formally rewrite the 2-point 
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function as 

W( Xl ,X 2 ) = 0(0, ^(T^e" 5 ^ 

= ^ Jv^vfoM^M^MHMMKl^^T}. (4) 

For the sake of the following discussion, let us stress that this gives the 2-point function 
in vacuum, (O|0(xi)0(a;2) |0). The same formalism can be applied to more general cases, 
such as for instance the 2-point function at finite temperature r, in which case instead of 
■00 [0] m (4) one should use the appropriate thermal state Vv[0]- A case of interest in the 
following is when the state is a coherent state of the theory. 

These are all well-known results from QFT that crucially use the presence of the 
Minkowski background. In particular, they can be extended to linearized GR around 
Minkowski, using the temporal gauge [19]. On the other hand, the situation is more 
subtle in full GR, where all the physical states have to solve the Hamiltonian constraint 
and thus have zero energy. Therefore the kernel can not have a decomposition like (2), 
but rather of the following type, 

K[gi,g 2 } = ^2*Pn[gi}^ n [g2}- (5) 

n 

Here n labels a complete basis of physical eigenstates, and it could be a continuous label (in 
which case the summation above is truly an integration). In this expression, the physical 
time is contained in the metric. 1 Consequently, the evolution equation (1) becomes a 
stability condition that any physical state has to satisfy, 

VW#2] = j V 9i K[g 1 ,g 2 ] ^ ph [#i]- (6) 

Equations like (5) or (6) encode the core of GR, and a successful theory of quantum 
gravity should be able to realize them explicitly. The spinfoam formalism provides us 
with a definition of the non-perturbative path integral which can be used to investigate 
(5) and (6). Here we provide a simple toy model in which we can realize them exactly. 
Furthermore, this explicit model will allow us to address another key issue, which concerns 
the interpretation of (4) in quantum gravity. As thouroughly discussed in [4], in full GR 
there is no special vacuum state singled out from the kernel, thus we can not apply the 
definition (4) straighforwardly. It is more sensible to consider the 2-point function on a 
semiclassical state ^[g], defined as a physical state peaked around a classical metric q. 
Then one can use a formula like (4), where instead of the ambiguous vacuum state ipo[g], 
one uses a semiclassical state, 

W q (x 1 ,x 2 ) = ^ J Vg 1 Vg 2 g 1 (x 1 )g 2 (x 2 )^ q [g 1 ]^ g [g 2 }K[g 1 ,g 2 }. (7) 



1 Assuming that (5) is peaked on a classical trajectory g^ between gi and <?2, the physical time can 
be computed evaluating the goo component along a geodesic, see discussion in [19]. 
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Upon suitable gauge-fixing, 2 this is a well-defined prescription, and indeed one can con- 
sider the 2-point function on a semiclassical state also in conventional QFT. However, this 
prescription requires the explicit knowledge of To be a physical semiclassical state, 
it has to (i) satisfy (6), and (ii) be peaked around a classical geometry q with minimal 
uncertainty. As discussed in [2], such a state should be peaked around both the intrinsic 
Rij(q) and extrinsic JCij(q) curvatures of q, so formally 

(Rij) = J Vg%[g]R tJ * q \g] = R tJ (q), (JC tJ ) = J Vg%^\t l3 * q [g] = JC tJ (q). (8) 

The latter are in general complicate equations to define and attempt to solve. To by- 
pass this difficulty, an interesting alternative is to find the analog of (3) in background 
independent quantum gravity, so to be able to extract the boundary state directly from 
the kernel. 

In the rest of this paper, we investigate these ideas in the simple toy model provided 
by a 3d Riemannian quantum tetrahedron. Using the spinfoam formalism, we give an 
explicit realization of (5) or (6), which we use to study the relation between the kernel 
and the boundary state, as well as its consequences on (7). 



3 The tetrahedron 

We consider the quantum tetrahedron introduced in 
[20, 3], see Fig.l. Given a tetrahedron of edge lengths 
£ e , we fix four opposite ones to j, and the remaining 
two to a and b. We orient it as in Fig.l in such a way 
that we can think of it as representing the evolution 
of the edge a into the edge b, in a "time" j. 




Figure 1: The tetrahedron. 



The dynamics of this model was studied at both 
the classical and quantum level in [20]. The classical dynamics is encoded in the Regge 
action S^(a,b,j) = ^2 e £ e 6e(£), where 9 e are the dihedral angles of the tetrahedron. The 
quantum dynamics can be studied a la Ponzano- Regge [14], making the ansatz that the 
lengths in the quantum theory can only have half- integer values £ e = j e + \ with j e G N/2, 
and associating with the tetrahedron the amplitude 



K[a,b,j] = y/d a d b 



a 
b 



(9) 



In this expression d a = 2a + 1 is the dimension of the spin-a representation of SU(2) and 
the {6j} is Wigner's 6j-symbol for the recoupling theory of SU(2). The simple form (9) 
allows us to study explicitly if a formula like (5) is indeed realized, and what its connection 
with semiclassical states might be. In particular, the stability condition for physical states 
(6) simply reads 

^&)=^X[a,MV(a)- (10) 



2 The issue of gauge-hxing is explicitly under control in the linearized theory [19]. There one evaluates 
the kernel and boundary states in the temporal gauge. Then the presence in (7) of the field insertions 
requires an additional time-independent gauge-fixing, such as the Coulomb gauge. 
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Before proceeding, notice that if we view this tetrahedron as part of a triangulation 
of flat 3d space between the two planes, we can also introduce the "asymptotic time" 

T(n\ + + 6V(a,b,j) 

nj) = t2-h= ab , 

where 

V(a,b,j) = ^-^Aa 2 b 2 j 2 - a 2 b 4 - a 4 b 2 (11) 

is the volume of the tetrahedron. In the isosceles case a = b = jo, we have T = 
a/ (2j 2 — j'o)/2 and V = j'qT/6. The latter quantities will frequently enter our calcua- 
tions below. Notice that fixing the bulk edges (to j) we did not fix the asymptotic time 
T, but we indeed fixed the way T depends on the (intrinsic and extrinsic) boundary ge- 
ometry. This fixing can then be viewed as the analog of the temporal gauge fixing used 
to define (5) in full GR, see [19]. 



3.1 Diagonalizing the kernel 

Consider the kernel as a dj-bj-dj matrix -Kab[j] — K[a,b,j]. Unexpectedly, this matrix 
satisfies K 2 = 1 for all j, thus the evolution generated by it is unitary. 3 As K 2 = 1, K is 
diagonalizable and all its eigenvalues are e n = ±1. Using the notation ip n (a) = (a\n) to 
indicate the n-th eigenvector in the basis a, we have 

2j 

K ab [j] = Y,£nMa)Mb)- (12) 

n=0 

This is an intriguing formula: in spite of the simplicity of the model, gives a non-trivial 
and explicit realization of (5) in quantum gravity. With respect to (5), notice the presence 
here of two possible eigenvalues ±1. This is the usual manifestation of the fact that the 
PR model sums over both orientation of the tetrahedron. 

The next step is to work out the eigenvectors solving the associate eigenvalue problem, 

J2K ab m(b) = ±^(a) (13) 

6=0 

It would be important to find the general solution to this problem, but we have only 
partially succeeded in doing so. It is easy to show that for any value of j, ipo(a) = \fd a is 
an eigenvector with eigenvalue e = 1. This can be done using the integral representation 
of the {6j}, 

J dgi dg 2 Xa(gi) Xb{92) Xjiom) XA9W2 1 )- ( 14 ) 

3 This might generate some confusion, as it is often stated that evolution in quantum gravity can not 
be unitary, due to the absence of a physical clock variable. However, the model here considered has a 
preferred clock variable built in it, and the evolution is unitary with respect to it. 



/ a j j \ 
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Then, using the fact that J2b Xb(g) — $(9), one can immediately show that J2 b W a bipo(b) = 
tpo(a). Similarly one can show that ipi(a) = (— \) a j ^[d~ a is an eigenvector with eigenvalue 
€1 = (— l) 2j . The remaining eigenvectors have to be found on a case by case basis for 
fixed j. For instance for the case j — 1 we have 

ipo(a) = (l,V3,VE), ipi(a) = (l,--J=,-^=) 

with eigenvalue 1, and ^2(0) = ^abc^o(b)ipi(c) with eigenvalue —1. 

For the general case, we can rewrite the eigenvalue problem (13) as an integral equation 
using the expression (14) as shown in Appendix A. However at this stage we do not 
have explicit solutions, so we can not study exact physical semiclassical states as linear 
combinations of eigenstates. In the next section we show how to solve it approximately. 



3.2 Semiclassical states 

The perturbative expansion we consider is the large spin limit (in which the half-integers 
approximate real lengths), where the {6j} symbol is dominated by exponentials of the 
Regge action, 

^ /10 t// u ■< cos ( g a[ a 'M + 7) • ( 15 ) 

y/12irV(a,b,j) v 4/ 

This is indeed the property that makes it possible to show that the quantum theory based 
upon the {6j} has a sensible semiclassical limit [14, 22, 3, 5, 23]. 

The lowest order of ty q should satisfy (13) with the linearized kernel and be peaked 
around q. In the simple setting considered here, the intrinsic and extrinsic geometry of q 
are specified giving a value j for edge length and a value 9 for its dihedral angle. We 
choose the latter in such a way that the complete background tetrahedron is given by the 
isosceles configuration a = b = jo- For the exterior dihedral angles associated to jo an d 
j, elementary geometry gives respectively 



cosg = - 4 / 2 2 3 j , cos0 = - ./° , 2 . (16) 



By analogy with the continuum, we expect (8) to be implemented in this approximation 
by a Gaussian state around q = (j , ), for which we make the following ansatz, 

%(a) = N exp{-^(a - j ) 2 } cos (^9 + </>). (17) 

Here iV is the normalization, a G R+ the width and is a phase that we leave un- 
determined for the moment. For a scaling linearly with the spins, this Gaussian is 
peaked on q in the large spin limit (e.g. [16]). This kind of states have been exten- 
sively used as ansatz for leading order semiclassical states in the recent spinfoam literature 
[2, 3, 4, 5, 6, 7, 13, 16, 21]. Yet, a semiclassical state typically has a precise width uniquely 
fixed by the dynamics (e.g. the factor u in the exponent of the vacuum Gaussian of the 
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harmonic oscillator). Thus a crucial question is whether also in spinfoams the width can 
be fixed requiring the state to solve the dynamics. Here we address the question in the 
specific model (13), and find a positive answer: as we show below, (17) is an eigenstate 
of the kernel for a unique choice of the width, given by 

4j 2 - jg = _jl cos9 > Q (18) 
8j Sj 



^2K[a,b,T\%(a) ~ #,(&) for = dfi + -(tt - 9), (19) 



and for two choices of 0, corresponding to eigenvalues ±1: 

.-0 + 

4 

a 

^2K[a,b,T\%(a) ~ -*,(&) for = + 1(tt - 0) + tt. (20) 

a 

Here ~ means at first order in the large spin limit. The presence of both eigenvalues 
could be anticipated from (12), and show that both approximate eigensolutions can be 
constructed as linear superposition of exact eigenvectors restricted to the same eigenvalue. 

Proving these results is straighforward, as in the large spin limit the eigenvalue problem 
(13), using (15) and the ansatz (17), amounts to a simple evaluation of Gaussian integrals. 
However, there is a clear restriction that should be kept in mind: formula (15) holds only 
if the volume V appearing in it is real, namely if V 2 defined in (11) is positive. For a 
generic configuration (a,b,j) this is not always the case. The quantum range of a and b 
is [0, 2j], and the condition V 2 > is violated when the endpoints are approached by one 
of the variables. Specifically, we have V 2 ~ for a and/or b close to zero, and V 2 < for 
a and/or b close to 2j. In both cases the asymptotics is not (15) (we report the relevant 
formulas in Appendix B for completeness). The calculations show below thus only apply 
in the regime a ~ b ~ j where we can use (15). 4 

We now follow a standard procedure (e.g. [16]), but with particular care about the 
phase. In the large spin limit, we can approximate the summation in (13) with an integral, 
and using (15) we have 



exp jie(s R (a, b,j) + ^ - -^(a - j ) 2 + ^(y#o + ^ 

Next, we expand the Regge action around the background q to quadratic order, obtaining 
S R (a, b, T) = 1d 3 Q + + ± 9 + - G lm 5j t 5 Jm , (22) 



4 A more thorough approach could also be possible: in [24] a unique expression that interpolates 
between the three regimes (V 2 > 0, ~ 0, < 0) was given. However it leads to a more complicate analysis, 
which we do not attempt here. 
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where we introduced the shorthand notation Sji = (a — jo, b — jo). The second derivatives 
of the Regge action where computed for instance in [5], and evaluated on q give Gu = 
G 22 = cos9/T and G 12 = G 2 i = -1/T, where we recall T = v / (2j 2 - j 2 )/2. 

The first three terms of (22) contribute to the overall phase of (21), which reads 



(23) 



The presence of a phase term in a exponentially suppresses the summation in (21). There- 
fore the configurations 77 = — e for which the a-phase vanishes dominate (21), and we can 
write 



AT f 

Y,K\a,b,T]m q {a) ~ - £ da 



d a db 



exp < le 



12ttV (a,b,j) 



x 



(24) 



§<?o + 2d,-6> + J 



4~( a - -^o) 2 + i| G lm 8ji 5j n 



This is a Gaussian integral whose leading order is obtained at 



d a db 



12nV(a,b,j) 



a=b=]o 



2 

ttT' 



(25) 



Then 



<l<i <>X P <! - -^(a - Jo) 2 + ^G /m 5j ; 5j r , 



(77T 



1 - 2ieaG 



11 



exp i 



where a' a complex function of a, e, Gu and G\ 2 that can be easily computed. In 
particular, requiring a' to be real implies Aa 2 = Aa' 2 = {G\ 2 — Gf^ 1 = T 2 /sm 2 9, 
independently of e. Plugging the values for the second derivatives of the Regge action 
given above we obtain the result (18). Using this result, a little algebra gives 



an 



v^ ex p{Ki~D}- 



1 - 2ieaG u 

The modulus of this quantity exactly cancels the contribution (25), so (24) reads 



(26) 



N 



J2 K fo b, T] *,(a) ~ — J2 ex P \ ie 



n-9 



-f9 + 2d J 9 + 9 



-^-j ) 2 \. (27) 



Thanks to the independence of a from e the summation in (27) gives a cosine and we 
finally get 

K[a, b, T] ,(o) ~ N exp{-^(6 - Jo ) 2 } cos (^9 + 2d 3 9 + ^ - </>) . (28) 



At this point it is easy to fix to obtain a result proportional to (17): the choice <fi = 
djO + {jx — 9)/ 4 leads to an eigenstate with eigenvalue +1, whereas the choice = djO + 
(7r — 9)/ 4 + 7r leads to an eigenstate with eigenvalue —1. With these choices of <fi and with 
a given in (18), (17) are approximate eigensolutions of (13) in the large spin limit. 

These results can be confirmed numerically, see Fig. 2. We consider the simple equi- 
lateral configuration j = j, for which notice that cos#o = cos9 = —1/3. 
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Figure 2: Numerical study of the ratio between (28) (with <f> giving the positive eigenvalue) and (17) 
in the equilateral jo = j case, with j = 80. On the left panel, a wide range of admissable values of a 
is shown. The behaviour is constantly around 1 only in a region with a around 80. For very small a, 
we enter the V 2 ~ regime, and (28) is not anymore equal to (17), as shown by the presence of the 
oscillations. Similarly the a ~ 2j regime shows a (still oscillating) exponential decaying. The region 
where the approximation is valid is shown in the right panel. 

Before concluding this section, some comments are in order. 

1. We started with a Gaussian with a generic width a. Requiring the stability of this 
Gaussian under the quantum evolution leads to a unique choice for the width. So 
the width can be fixed dynamically. This is independent of the phase term in (17). 

Furtermore, notice that this choice is precisely the one made in [3, 5] in order to 
obtain the correlation of an harmonic oscillator. Here we show that this choice arises 
naturally from a dynamical requirement. 

2. The choice of a real phase (the cosine) is rather natural given the form of the kernel 
asymptotics (15). On the other hand, a complex phase of the type exp{i9dj/2} 
has been more commonly used in the literature. This choice has the advantage 
of selecting only one of the two exponentials in (15), thus introducing a notion of 
"observer-induced" orientation. 

Concerning this choice, the reader can easily check that an ansatz of the type 

*,(o) = N exp{--^(a - j f + i^-9 + i<f>} (29) 
4a Z 
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can not solve the eigenvalue problem (13). On the other hand, one can show pro- 
ceeding as above that (29) with a given in (18) and in (19) satisfies the equation 



Y,K[a,b,T\* q {a) = * q {b). (30) 

a 

We can read the above equation as a stability condition where the state is turned into 
its complex conjugate. This can be also compared with the semiclassical tetrahedron 
used in [16]. 

3. As shown in [5], the modified Bessel functions of the first kind I n (z) satisfy 

y/I Q {Aaj Q ) - / 2jo +i(4(7jo) 4a 

This means that also this Bessel state approximately solves the eigenvalue problem. 
However, it does not solve it exactly. The advantage of this boundary state is that 
it admits a simple SU(2) Fourier transform, which was used in [6] to cast the 4d 
graviton into an integral representation. 

The consequences of choosing the physical boundary state computed here in the eval- 
uation of the graviton propagator are discussed below in Section 5. 



4 Physical state in the general boundary formalism 

In this section we discuss a more speculative construction. Above we computed the leading 
order approximation of a state that (i) solves (6) and (ii) is peaked around a semiclassi- 
cal geometry, but we were not able to extend this result to the non-perturbative theory. 
However, the situation becomes somewhat simpler in the so-called general boundary for- 
malism, where also the four bulk edges are varied freely. Then the boundary state has 
to carry information on the background value of the (intrinsic and extrinsic) geometry of 
all six edges, and correlations between all six of them can be computed. Such a general 
boundary formalism has been advocated to be the most natural area for quantum gravity 
[25], and is indeed the setting used for the 4d spinfoam graviton calculations [2, 4, 6, 8]. 

In general we might expect the general boundary problem to be more complicated 
to deal with than the fixed-time setting. Yet in the particular 3d case the fact that the 
theory is topological strongly simplifies the analysis, because once the topology and the 
triangulation are fixed, there is a single physical state. Assuming trivial topology, the 
latter is given in the group representation by 

* = n ^ (32) 

ind/ 

where gj = Yliedfdi represents the gravitational holonomy on a closed path df. The 
product is over the independent faces, and the condition F = is ensured everywhere. 
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The spin representation of (32) can be obtained projecting it on the spin network basis 

* {jl} (gi,...g 6 ) = l[inl[D^(g l ). (33) 

n I 

In the case considered here, the triangulation is given 
by a single tetrahedron, and (32) is a product of four deltas 
ensuring the flatness of each face. Only three faces are 
independent, so we can get rid of one delta. We fix the 
orientation of the edges as in the adjacent figure. Choosing 
to neglect ^(gig^gs) we can simply write the state 




1 

M9i, ■ ■ ■ 9e) = 5{g,g 2 gS) 5(g 5 g,g 6 ) S^g^ 1 ). (34) Fi S ure 3: Orientation. 



This state is physical: it satisfies F = everywhere and it is thus a solution of the 
Hamiltonian constraint. 

To obtain the spin representation of this state, we project it on the spin network basis 

*{j,}(0i,---06) = Yl D minM ■ ■ ■ D mjn 6 (96) (mim 2 n 3 ) (nim 5 n 6 ) (men 4 n 2 ) {m A n 5 m 3 ), 

rrii,ni 

where the (m 1 m 2 m 3 ) are normalized Clebsch-Gordan coefficients. A straightforward cal- 
culation then gives the projection 

MJi ■ ■ ■ k) = d 9l ... dg 6 $ (g u . . . g 6 ) # 0i} (<?i, . . . g 6 ) = \ {} f f 1 . (35) 

J [ J4 j5 J6 J 

We see that the physical boundary state coincides with the kernel, so that correlations 
(see next Section) now read: 

This result might look surprising at first, but it can be understood as follows. The 
boundary state V'oO'i • • -3%) is the state induced by the (exterior) bulk geometry onto the 
tetrahedron. We have assumed a trivial boundary topology (homomorphic to S 2 ) and 
a trivial bulk topology. Therefore we obtained a spinfoam amplitude in {Qj} 2 which is 
naturally associated to the triangulation of the closed S 3 manifold with two tetrahedra. 
In this context it is a natural result. 

Notice that to go back to the time-gauge setting used in the previous Section it suffices 
to gauge-fix four opposite representations of the spin network state to j. We can interpret 
this gauge-fixed spin network state with the four representations fixed to j as if it was 
induced by a non-trivial bulk topology. Indeed, we can imagine topological defects on 
each of these edges inducing non-trivial holonomies around these edges. These are in 
turn interpreted as particles traveling along these edges and j is the proper time along 
the particle trajectory [26, 22]. From this perspective, the gauge-fixed spin network state 
with fixed j appears as a physical state for a non-trivial topology, or, equivalently, for the 
coupled system "gravity + particles". 
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5 Graviton propagator revisited 



The results presented so far give new inputs to the physical interpretation of the 3d 
graviton propagator computed in [3, 5]. Let us recall that in the time-gauge setting, (7) 
is realized as [3] 

W(joJ) = \^J2Hja)Hj b )^ q [a,b]K ab [j] (37) 



a,b 

where h(J a ) = j a (ja + 1) — Jo (jo + 1) represents a field insertion. This quantity gives 
the quantum correlator between the two opposite edges, and formally can be considered 
as a projected component of the graviton propagator. The kernel is (9), 5 and in [3] the 
generic Gaussian ansatz with a complex phase (29) was used. Here we want to study what 
happens when we use the physically determined unique boundary state. The leading order 
in the large spin limit can be easily computed following the same procedure of Section 
3, and noticing that h(j a ) ~ 2j 5j a - The result is a phased Gaussian integral whose 
evaluation gives 

w( . 4 E £=± exp{*e[f + 2^-20]}(4- 1 ) 12 / v ^O: 

U0,3) fo £ e=± exp{ie[f + 2d jo 0-20]}/v/detX ~ 1 J 

where the quadratic matrix is 

A e = — [ -i n sm \ Q , det A e = - - . - e te{0 a). (39) 

Simple algebra shows that exactly cancels the phases in (38) for both choices (19) and 
(20), and finally we get 

W( 3o ,j) = — co S 9. (40) 
Jo 

Following [3], this can be interpreted as the real part of the correlator of an harmonic 
oscillator of frequency u = 9/T. However, the fact that we are now working with the 
real part allows us to have a direct interpretation in terms of the inverse distance traveled 
along the bulk edge. Indeed, using the explicit values (18) and (16) for a and cos# we 
have 

W(joJ) = -±:, (41) 

independently of jo (although the very definition of the leading order approximation is 
based on the presence of the background jo). This remarkable result shows explicitly at 
leading order that the use of a physical boundary state leads to the correct behaviour of 
the correlations. 

Yet to make a more precise connection between this correlation and the (projections 
of the free) graviton propagator, there are two crucial issues that have to be addressed: in 
the continuum theory the 3d propagator is (i) gauge-dependent, and (ii) is a pure gauge 



3 Up to factors cL irrelevant for the leading order analysis. 
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quantity, consistently with the absence of local degrees of freedom of 3d GR. These issues 
were not addressed in [3, 5], and the question remained open for some time. 

The answer to point (i) is rather simple: the expression (37) is defined in the time- 
gauge obtained fixing the bulk edges to j. Consistently also the physical boundary state is 
evaluated in the time gauge. Now in the continuum the time gauge is only a partial gauge 
fixing (e.g. [19]), and one needs a further spacelike gauge fixing (such as the Coulomb 
gauge) in order to compute the the 2-point correlation, which involves the insertion of the 
gauge-dependent quantities h a b(x). In (37) on the other hand, the correlations are studied 
via the gauge- independent observables lh (they are SU(2) Casimirs), thus the time gauge 
is a full gauge-fixing, and the correlation (41) depends only on it. 

Answering point (ii) is on the other hand more subtle. This question was addressed 
in [9], where it was shown in 3d Regge calculus on a infinite regular lattice that the 
propagator is a pure gauge, namely that the propagating quantity in the correlations 
comes entirely from the gauge-fixing term in the discrete action. It would be useful to 
extend those results to the spinfoam formalism considered here, however that calculation 
relies on the fact that one is working on an infinite lattice, and thus gauge degrees of 
freedom can be identified as translations on the lattice. This is not possible from the 
viewpoint of a single tetrahedron. As shown in [9], the left over of gauge invariance on 
a single tetrahedron can be identified at the perturbative level as the global symmetry 
corresponding to the following global rescaling on the background £ e , 



When we choose the time gauge we clearly break this symmetry. However one has to show 
that it is precisely this gauge-fixing term that propagates in (41), and nothing else. To 
do so, we consider the general boundary formalism discussed above, where the symmetry 
(42) is present. 

In the general boundary formalism, using the physical boundary state (35) leads to 
the formula (36) for the correlations between edges. The perturbative expansion of this 
quantity can be studied again using (15), which gives in the large spin limit 



The constant first term in the square bracket shows immediately that the summand is 
divergent, consistently with the fact that we have not fixed the gauge. Furthermore, 
notice that the dynamics is only encoded in the oscillatory second term. Imagine now to 
gauge-fix the symmetry (42) introducing a la Faddeev-Popov terms like exp{— 5j 2 }. If this 
gauge-fixing term has a non-diagonal structure then it will give the dominant contribution 
to the correlations, with the Regge action in (43) exponentially suppressed. In this case, 
we can clearly claim that at least perturbatively the correlations are pure gauge. If on 
the other hand the gauge-fixing term has no off-diagonal structure, then the correlations 
would indeed come from the dynamical Regge term, but because of its oscillatory nature 
they would be completely suppressed at large scales. They would then be some sort of 
quantum noise of no semiclassical counterpart. 



5£ e i-> 5£ e + N£ e . 



(42) 




(43) 
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6 Conclusions 



We considered the stability condition for physical states in quantum gravity. Within the 
simple toy model of a 3d Riemannian tetrahedron, we realized explicitly this condition as 
an eigenvalue problem in the spinfoam formalism, given in (13). We wrote this problem 
as an integral equation (see Appendix A), and found some exact solutions. At leading 
order in the large spin limit we were able to find a class of solutions which can also be 
interpreted as semiclassical states for the linearized theory around a fixed background. 
This class of solutions are Gaussian wavefunctions (17) with a width that is dynamically 
fixed to the value given in (18) by the stability condition. This result positively answers 
the question that one can fix the width of this semiclassical states looking at the spinfoam 
dynamics. Furthermore, in this specific toy model, the amplitude fixed in this way turns 
out to be exactly the one giving correlations which are at leading order in agreement 
with the discretization of the continuum linearized quantum theory. Finally, our results 
show that the PR model supports the dynamical propagation of stable wave packets of 
geometry. 
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A Analyzing the {6j} kernel 

We start from the integral expression (14). Expressing the group elements g± t 2 in term of 
their class angles 0i )2 and their rotation axis G S 2 , this becomes: 

J J'} = ^/ #i#2ffln 2 0ism 2 02Xo(0i)X6(02) J ^dyXj(<t>-)Xj(<l>+), (44) 
with the mixed class angles defined in term of y = Ui ■ u-i as: 

cos0-|- = cos 0i cos 02 =F 1/ sin 0! sin 02- (45) 
With y = cos a, notice this can re- written as: 

cos0_ = cos 2 — cos(0i — 2 ) + sin 2 — cos(0i + 2 ), 
cos0 + = sin 2 — cos(0i — 02) + cos 2 — cos(0i + 02), 

which shows that the angles <fi± run between (0i — 2 ) and (0 X + 2 ). Notice also the 
constraint between the angles: 

COS 0_ + COS + = cos(0i — 2 ) + cos(0i + 2 ) = 2cOS0iCOS0 2 . (46) 
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Using the relation between the SU(2) characters and the Chebyshev polynomials of the 
second type, Xj{4>) = sin(2j + 1)0/ sin0 = [/jy(cos0), we can re-express the {6j} as an 
integral of products of Chebyshev polynomials over x^ 2 = cos0i 5 2 and y. For instance, 
we have: 

| ° \ j | = ^y dx 1 dx 2 dx-dx + 5(x- + x + - 2x 1 X 2 )U2a(xi)U2b(x2)U2j(X-)U 2 j(x + ). 

Expression (44) suggestes to study the eigenvalue problem (13) in terms of the following 
Fourier transform of ip a , 

^(0) = ^V / ^Xa(0)^a. (47) 

The eigenvalue equation then becomes a convolution equation, 



V>(0) = ±/ #2 sin 2 cf>2 X(0,0 2 ,j)V(0 2 ), (48) 

JO 

with kernel given by 

K{4>,<h,j)= [ #i sin 2 faKfa, 2 )?(0,0i), (49) 
./o 



where 



^(01,02) = J dyXj(<f>-)Xj(<f>+), <K0,0l) = ^2d a Xa(<f>i)Xa(<f>)- 



(50) 



K(4>i,(p 2 ) can be evaluated exactly writing it in term of Chebyshev polynomials while 
5(0, 0i) is a distribution over SU(2) which can be directly evaluated using the explicit 
expression for the characters, Xa{4>) — sin(2a + 1)0/ sin0. 

Exact solutions to the convolution equation are still under investigation. At least, we 
already know that gaussian states on SU(2) will approximatively solve it in the asymptotic 
limit for large spins since they are the Fourier transforms of the gaussian states on the 
representation labels. 

B Extremal cases 

The first and last rows of the matrix if a &L?] admit simple expressions, that we report here 
for completeness and future use. Notice that these are cases for which V 2 is not positive 
so the asymptotics (15) do not hold. 

Case 6 = 



In this case we have V 2 = and 



j J J ~" d 3 



(51) 
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Therefore gj(a) = -Kao[j] = (— l) 2 - J+a ^jS. This is a rather simple expression, plotted on 
the left panel of Fig. 4. 

For generic V 2 ~ configurations, the asymptotics of the {Qj} is [14] 



\/64S 



cos $ M(z) + sin $ Bi(z) 



(52) 



where S is the four times the product of the four triangle areas, Ai and Bi are Airy 
functions [27], and 



-(3V) 2 S~I } V 2 >0 
(3\V\) 2 S~I, V 2 <0 



(53) 



Case b = 2j 

In this case we have V 2 



< o and 



a J J 
2j j j 



Aj 2 T{2jf 



p(a,j)T(2j + a)T(2j-a) 



(54) 



where 



p(a,j) = (2j + a){2j - a){2j + a + 1). (55) 

The function fj(a) = K a2 j[j] has a more interesting behaviour, and a peak approximately 
at a/j/2, see the right panel of Fig.4. The position of the peak can be computed analyt- 
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Figure 4: Left panel. A typical plot for gj(a). Here j = 50, and for clarity only part of the 
full range a G [0,100] is shown. Right panel. A typical plot for fj(a). Here again j = 50 and 
a e [0, 100]. The peak is at = 5. 

ically: 

d aWj (a) = -Aj 2 y4jTTr( 2j ) 2 «M)-& + l)P(">MM23-a)-M23 + a)] (5g) 

V2a + lp(a,j) 

where 

q(a,j) = 5a 3 + 8j 3 + 6a 2 (j + 1) + 2a - Aaj(j - 1) (57) 
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Mn) = -lE + J2l (58) 

k=l 

is the digamma function. Its power series for large values of the argument starts with 
ln(n — 1) + l/2(n — 1). Using this approximation, the leading order for big j of the 
numerator of (56) vanishes at 

2j - a(4a + 3) = (59) 

from which 

U ^ /J 

a = 

r2 



a =-^ V 9 + 32j-3j~W|. (60) 



For generic V < configurations, the asymptotics of the {6j} is [14] 

m ~^mv cosie ^ (61) 

with 

* = 5>'e-^)Re0e, fi = ^J e Im^ e . (62) 
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